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ABSTRACT 

In this paper, the third of a series dedicated to the investigation of the nuclear 
properties of spiral galaxies, we have (i) modelled the WFPC2 F606W nuclear surface 
brightness profiles of 41 spiral galaxies presented in Carollo et al. 1997c, 1998 with 
the analytical law introduced by Lauer et al. 1995, and (%%) deprojected these surface 
brightness profiles and their analytical fits, so as to estimate the nuclear stellar 
densities of bulges of spiral galaxies. We find that the nuclear stellar cusps (quantified 
by the average logarithmic slope of the surface brightness profiles within 0.1"-0.5") are 
significantly different for i? 1 / 4 -law and exponential bulges. The former have nuclear 
properties similar to those of early-type galaxies, i.e. similar values of nuclear cusps 
for comparable luminosities, and increasingly steeper stellar cusps with decreasing 
luminosity. By contrast, exponential bulges have (underlying the light contribution 
from photometrically distinct, central compact sources) comparative shallower stellar 
cusps, and likely lower nuclear densities, than i? 1 / 4 -law bulges. 

Subject headings: Galaxies: Spirals — Galaxies: Structure — Galaxies: Fundamental 
Parameters — Galaxies: Nuclei 
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1. Introduction 

The galactic nuclei are the repositories of low angular momentum material sunk to the 
centers over the lifetime of the parent systems. Therefore, they are likely to hold answers to 
important questions related with the origin of structure in the parent galaxies. In this perspective, 
establishing the demographics of galactic nuclei along the entire Hubble sequence lies at the heart 
of our understanding of the complex process of galaxy formation and evolution. 

Observations of nearby ellipticals and lenticulars with the Faint Object Camera (FOC), the 
Wide Field Planetary Camera (WF/PC) and the Wide Field Planetary Camera-2 (WFPC2) aboard 
the Hubble Space Telescope (HST) have revealed that the nuclei of these galaxies are complex 
environments (e.g., Crane et al. 1993; Jaffe et al. 1994; Lauer et al. 1995, hereafter L95; Forbes et 
al. 1995; Carollo et al. 1997a, hereafter C97a; Carollo et al. 1997b, hereafter C97b; Faber et al. 
1997, hereafter F97). They show surface brightness profiles that increase down to the innermost 
point measurable at HST resolution, i.e., I(r) oc r~ 7 as r — > (where I[r) is the surface brightness 
at the radius r, and 7 > 0); furthermore, several galaxies host stellar and gaseous disks, unresolved 
nuclear spikes, double nuclei. These inner features might possibly be related to the presence of 
massive black holes (e.g., Lauer et al. 1996). 

By contrast, much is still to be learned about the nuclear properties of nearby spiral galaxies 
at HST resolution. F97 found that the surface brightness profiles of the three Sa-Sb bulges 
present in their sample show a behaviour similar to that of early-type spheroidals of comparable 
luminosity. The same result was found by Phillips et al. (1996) for the three spirals of type earlier 
than Sc contained in their WF/PC F555W sample of 20 disk galaxies. Furthermore, Phillips et al. 
found that later type spirals show instead (almost) flat nuclear profiles, and suggested that the 
nuclear properties of disk galaxies are more closely related to those of nucleated dwarf galaxies 
than to those of elliptical galaxies. Further exploration is necessary to assess how the nuclear 
properties scale with the properties of the spheroidal component. This is likely to provide feedback 
information about the epoch and processes of nucleus, bulge, and, ultimately, galaxy formation. 
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In order to address these issues, we have performed a WFPC2 snapshot survey in the F606W 
filter of the nuclei of 107 (mostly Sa to Sc) disk galaxies. In paper I (Carollo et al. 1997c) and 
paper II (Carollo et al. 1998) we have presented the 75 targets imaged so far within our program. 
Our analysis shows that bulge-like structures are present in most of the galaxies. While in some 
cases these are "classical", smooth, featureless i? x / 4 -law bulges, in others they are better fitted 
by an exponential profile (see also Courteau, de Jong & Broeils 1996, and references therein, for 
similar results). The exponential bulges include two classes of objects: (i) dwarf-looking systems, 
whose surface brightness profiles within ~ 15" are well fitted by a single exponential. These 
galaxies are strongly bulge-dominated; their surrounding, faint regions ("disk/halo") show no 
signs of spiral arms, and have typically a quiescent, i.e. non star forming, appearance, (ii) Small 
exponential bulges embedded in dominant (spiral- armed/star- forming) surrounding disks, i.e., the 
inner exponential structures of double-exponential fits to the surface brightness profiles within 
~ 15". The exponential bulges as-a-class are statistically fainter than the featurless, smooth i? 1 / 4 
bulges, for constant disk luminosity and Hubble type. Resolved, central compact sources are found 
in most of the exponential bulges; the hosts of central compact sources often contain a barred 
structure. The nature of these compact source, and in particular their relation with e.g., star 
clusters and Seyfert 2 nuclei, is discussed in Carollo (1998). 

In this paper, the third of the series, we investigate the relation between the nuclear structure 
of spiral galaxies and the physical properties of inner disks and/or bulges. In particular, we 
(i) present the results of modeling the nuclear surface brightness profiles with the analytical 
law introduced by L95 (for the 43 galaxies of paper I and II for which we could perform the 
measurements), (ii) deconvolve the surface brightness profiles and their analytical fits in order to 
estimate the nuclear stellar densities, (Hi) study the nuclear properties as a function of the global 
properties discussed in papers I and II (e.g., i? 1 / 4 against exponential bulges), and (iv) compare 
the nuclear properties of our sample with those observed in early-type galaxies. The paper is 
organized as follows. In section 2 we briefly summarize the properties of the sample, the data 
used in our investigation, the procedure adopted for the data reduction, and the steps performed 
to derive the surface brightness profiles. In section 3 we present the results of the analytical 
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fits applied to the surface brightness profiles, and of deconvolving data and models in spherical 
symmetry. In section 4 we investigate the dependence of the nuclear properties on global galactic 
properties. We conclude in section 5. 

2. The Sample, the Data and the Data Reduction 

We address to papers I and II for details on the selection of the total sample, the data, and 
the procedure adopted for their reduction. Below we briefly summarize the important points. 

We selected a total of 92 spiral galaxies (from the UGC catalog for the northern hemisphere, 
Nilson 1973, and from the ESOLV catalog for the southern hemisphere, Lauberts & Valentijn 
1989), which satisfy the following constraints: (i) Angular diameter larger than 1 arcmin; (ii) 
Regular types Sa, Sab, Sb, and Sbc; (Hi) Redshift < 2500 km/s, to guarantee a high angular 
resolution in physical size; (iv) An inclination angle, estimated from the apparent axial ratio, 
smaller than 75 degrees. Fifteen additional galaxies were chosen among E/SO and SO/Sa, so as to 
have a reference sample for comparison with the literature on early-type galaxies. In this paper 
we use the surface brightness profiles of the 18 objects presented in paper 1, and of the 25 objects 
presented in paper 2, for which the light distribution was regular enough to allow us to perform 
the measurements. Some parameters for the total sample of 43 galaxies studied in this paper are 
given in Table 1. We adopt H Q = 65 km/s/Mpc throughout this paper. 

For each galaxy we acquired two WFPC2 F606W exposures of 400 and 200 seconds, respectively. 
The observations were taken with (gain 15 and) the galaxy nucleus centered on the PC, which 
has a scale of 0.046" px _1 , and a useful field of about 35" x 35". The observations were carried 
out in fine lock. The data reduction included: (i) the reprocessing of the raw images with the 
standard WFPC2 pipeline in order to use the most recent reference frames for flat fielding, bias and 
dark current subtraction; (ii) the addition and cosmic-ray cleaning of the two images by using 
the IRAF/STSDAS task "crrej"; (Hi) the removal of remaining hot pixels by interpolation; (iv) sky 
subtraction by determining the sky values from the WFC images, in areas furthermost from the 
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nucleus. The derivation of the surface photometry was carried out in IRAF. For the sample of 
paper I, we used two different ellipse-fitting programs ( "Galphot" , Franx, Illingworth & Heckman 
1989, and "Ellipse", R. Jedrezjeski), and two different approaches for the dust- and star-masking, 
so as to have an estimate of the uncertainties introduced by the complex dust and star forming 
structures. Both approaches only correct for patchy dust absorption and not for any smooth and 
extended dust component, which remains undetectable in our data. For the sample of paper 
II, we only used the task "Ellipse", since the previous study shows that the two packages give 
equal results within the errors. Absolute photometric calibration was obtained by applying the 
zero-point for the F606W filter provided by Holtzman et al. (1995), i.e. 

mag F 606W = -2.5 \og(counts / sec) + Z WFPC ,F606W (1) 
+51og(0.046) +0.1 

with the zero-point Zwfpc,F606W equal to 22.084. The subsequent term accounts for the pixel 
size. The constant shift of 0.1 mag was added following Holtzmann et al. in order to correct for 
infinite aperture. This calibration provides a rough conversion to V- Johnson magnitudes (i.e., in 
average within s=s 0.2 magnitudes, given the difference in bandpass between the Johnson V and 
the F606W filters, and the lack of a color term for our sample). The calibrated surface brightness 
profiles were corrected for Galactic extinction by using the values published by Burstein & Heiles 
(1984). 

3. Analytical Fits to the Nuclear Surface Brightness Profiles 

One of the main purposes of any analytical description of the light detected from astronomical 
sources is to provide a simple way to quantify the most relevat features and their variations from 
object to object within the same family of sources, and to compare these features among objects 
of different families. For the early-type galaxies, a popular description of the surface brightness 
profile has been proposed by L95, and used by Byun et al. (1996) to describe the nuclear properties 



- 7- 



of kinematically normal early-type galaxies, and by C97a, C97b to describe those of elliptical 
galaxies with kinematically-distinct cores. It reads: 



/( r ) = 2^y-i b (^r \i + (^-H (7_/3)A \ (2) 



r 



The parameter i4 is the break radius at which, in giant elliptical galaxies, the profile flattens 
to a more shallow - but always non-zero - slope, 7 is the slope of the inner profile (i.e., I(r) ~ r~ 7 
as r — > 0), /? is the slope of the outer profile, a controls the sharpness of the transition between 
inner and outer profile, and is the surface brightness at R\>. The profile of equation (2) provides 
an adequate description of the features present in the surface brightness profiles of early-type 
galaxies along their entire luminosity sequence. 

The structure that we detect in the nuclei of several disk galaxies requires a clear definition of 
what are the relevant features in these galaxies that we aim at quantifying. In particular, for the 
exponential galaxies which host a perfectly centered compact source, there is a conceptual choice 
which has to be made, namely whether to include the central source in the analytical description 
of the galaxy (in a few galaxies, the "central" compact source reported in Table 1 is actually 
slightly offset from the center of the outer isophotes. Therefore, it has not been considered in the 
isophotal fitting, and does not appear in the surface brightness profile). In principle, these systems 
could be described by an exponential profile modified to include a steep nuclear cusp slope. In 
order to check whether this is a physically acceptable solution, we fitted the exponential-bulge 
galaxies hosting a central compact source with a light profile described by: 



I{R) = Io [I + -^J exp(-R/R,). (3) 

For radii R much smaller than the "cusp" radius R c , this profile describes a cusp with slope 7. Iq 
is the central brightness of the esponential component, and R s the esponential scalelength. No 
acceptable (i.e., physically meaningful) fit could be found for any galaxy except ESO 499 G37. 
This demonstrates that "steep-cusp esponentials" could indeed exist, and that our procedure 
would be able to identify them. However, the fact that they generally do not provide a good or 
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meaningful fit supports the interpretation, based on morphology, that central compact sources 
are photometrically distinct components from their surroundings, and not simply a steepening of 
the light profile. On this basis, we considered these sources as an additional contribution to the 
to-be-fitted, underlying galaxy light. The radial profile of the latter is what can be meaningfully 
compared to the radial profile measured for the early- type systems. Indeed, the same choice was 
made by Byun et al. (1996), who excluded from their fits the light excess from the resolved nuclei 
embedded in the bulges of a few galaxies of their sample (e.g., NGC 1331, NGC 3599, NGC 
4239). Therefore, in order to perform a direct comparison with elliptical and lenticular galaxies, 
we adopted equation (2) to describe the nuclear properties of spiral galaxies. Some details on the 
analytical "L95 fits" are given in Appendix A. 

There is an intrinsic uncertainty in disentangling the contribution of the nuclear compact 
source from the unknown underlying galaxy light. In order to quantify this uncertainty, we 
performed several fits to the same galaxy, by considering different inner radial cutoff in its surface 
brightness profiles. The resulting inward extrapolations of the best fits to the outer points may 
differ significantly, but are likely to bracket the true shape of the underlying galactic light, free 
from PSF-blurring and central compact source contamination. We were able to fit 41 of the 
43 galaxies (for the two remaining objects, ESO 205G7 and NGC 3177, no acceptable fit was 
obtained). In Figure 1 we show the surface brightness profiles, and the fits with the Point Spread 
Function (PSF) convolved models, for the 41 galaxies. In Table 2 we list, for each galaxy, radial 
range adopted for performing the fits, and the corresponding best fit values for Rb, a, (3, 7 
and Hb = —2.5 log /ft. For those objects for which changing the radial range for the fits lead to 
significantly different solutions, we show in the table and in the figure the two "bracketing" fits. 

For each galaxy, the analytical fits to the surface brightness profiles were used to derive the 
average logarithmic nuclear slopes (7} corrected for PSF and central compact source effects. The 
average slopes provide a global representation of the underlying projected galactic light within a 
defined radial range. This was taken equal to 0.l"-0.5" (similarly to what was done by L95 and 
C97a for the early- type systems), which corresponds to physical scales smaller than w 100 pc 
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at the average distance of the galaxies (« 25 Mpc). In contrast to the values of the individual 
parameters in the best fits, which for the same galaxy can differ significantly from one best fit to 
the other (obtained within different radial ranges), the (7} values provide a robust description 
of the nuclear surface brightness profiles contributed by the underlying galactic light. For the 
galaxies for which two fits are presented, we adopted the arithmetic mean of the corresponding 
values of (7} . An estimate of the uncertainty associated with these measurements is given by the 
difference between the two values, and is typically smaller than ~ 0.2. The values of the average 
logarithmic slopes inside 0.1"-0.5", (7), are also listed in Table 2. 



3.1. Deprojections 



The analytical fits to the surface brightness profiles were deprojected in spherical symmetry 
by means of an Abel inversion (Binney & Tremaine 1987), to derive the luminosity density 
profiles p mo d(r)- The same was done using directly the data points [pdata{f)\- Despite the fact 
that the deprojection amplifies the fluctuations and errors present in the data, and although 
the deprojection of the data introduces systematic errors due to the effects of the WFPC2 PSF, 
at radii of ~ 0.1" (where the effects of the WFPC2 PSF are modest) the comparison between 
Pmod(f) an d Pdata(r) provides a rough estimate for the effect of the central compact source on the 
underlying stellar density, and more generally an estimate for the uncertainty associated with the 
measurements. 

In order to quantify the general accuracy of this estimate, and in particular the effects of 
the WFPC2 PSF, we carried out the following test. For the galaxies which do not host a central 
compact source, we started from their (PSF-free) L95 models, deprojected both these "intrinsic" 
profiles and those obtained from them after convolution with the WFPC2 PSF, and compared the 
corresponding densities at 0.1". The deprojections of the intrinsic and PSF-convolved profiles 
yielded densities at 0.1" (or 10 pc) within 10-15%. For the exponential galaxies which host a 
central compact source, we added to their (PSF- and compact source- free) L95 models a central 
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source described by a Plummer model of varying radius, in order to also assess the effect of 
such a light concentration on the central densities. We then repeated the same steps described 
above, i.e., deprojected these intrinsic profiles and those obtained after PSF-convolution, and 
compared the resulting densities at 0.1". In this case, the comparisons showed that for the range 
of compact source luminosities and radii observed in our sample, the straight deprojection of the 
data, ignoring the effect of the HST PSF, can be in error by up to a factor 50% (see Figure 2, 
where we plot, as an example, the results for ES0482G17). This error is much smaller than the 
effect of the compact source itself, which typically increases po.i", pio pc of a factor ten or more. 

The deprojected luminosity density profiles are shown in Figure 3, where the deprojections of 
the fit (pmod) are plotted as solid lines, and those from the data (pdata) as dashed lines. When two 
analytical descriptions were available for the same galaxy, in Figure 3 we plot the deprojection 
corresponding to the first of the two fits in Table 2. The values of the stellar densities at 0.1" and 
at 10 pc derived from deprojecting the data (pdatafi.v and pdata,i0pc, respectively) and the models 
(Pmodfi.i" and Pmod,io P c respectively) are given in Table 3. 

In the reminder of this paper, we use the values of Pdatafi.v-, Pmodfi.i" and (7) to describe the 
nuclear properties of disk galaxies. The choice of using the densities at 0.1" rather than those at 
10 pc is justified by the fact that the latter, although more physically motivated, often probe a 
radius smaller than the typical FWHM of the PSF. Nonetheless, we have checked that our main 
results remain unchanged when the densities at 0.1" are substituded with those at 10 pc. For 
those objects for which two values of p mo dfi.i" were available, we adopted the smallest of the two. 
This choice maximizes the estimate for the associated errorbar, since in these objects typically 

Pdata ,0.1" > Pmodfi.i"- 

4. Discussion 

The HST exploration of the nuclei of nearby early-type galaxies has highlighted the 
diversity between low-luminosity, rotation-supported, radio-quiet, disky-isophote galaxies, and 
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high-luminosity, pressure-supported, radio-loud, boxy-isophote systems. The diversity holds in 
fact down to the tens-of-parsecs scales: the former have power-law (steep-cusp) light profiles and 
high stellar densities; the latter show a break radius inside which the profile flattens to a shallow 
(but non-zero) cusp slope, and have low-density nuclei (e.g., L95). It is unclear however whether 
the connection between small and large scale properties in early-type spheroidals is operated by 
nature, i.e. it is the outcome of a different formation epoch, time scale or physical mechanism, or is 
operated by nurture, i.e. it results from subsequent galactic evolution. In the orthodox framework 
where bulges of spiral galaxies are the extrapolation of the spheroidal family toward fainter 
magnitudes, and in the light of the variety of structural properties shown by the central spheroidal 
structures embedded in disk galaxies (e.g., Kormendy 1993; Courteau et al. 1996; papers I and II), 
the differences and similarities between the nuclear properties of (i) i? x / 4 -law versus exponential 
bulges, and (ii) bulges versus higher luminosity spheroidals, can provide important feedback 
toward understanding the processes involved in the formation of spheroidals and their nuclei. 

In Figure 4 we plot the average nuclear slope inside 0.l"-0.5", (7), versus the total F606W 
magnitude of the spheroidal component Mfsogw- Triangles represent galaxies fitted by a single 
exponential, pentagons those fitted by two exponentials, circles the i? 1//4 -law bulges, and asterisks 
the early-type galaxies of L95, F97 (small symbols), and C97a, C97b (large symbols). Only 33 of 
the 41 bulges for which we could derive the mesurements of (7) appear in the figure: for seven 
of the remaining eight objects, no measurement of Mpm&w was available (see Table 1, galaxies 
labelled by "NGF" in the last column - ESO205G7 is not included in these seven object, since for 
it neither Mfqqqw nor (7) could be derived); furthermore, we excluded NGC 488 from the plots 
since its measurement of Mp606W is uncertain (although it definitely shows an it^-law profile, 
see Figures 6 and 8) . We verified that the inclusion of NGC 488 in our analysis does not change 
our conclusions. 

There is a clear dichotomy between i? 1 / 4 -law and exponential bulges, which holds even when 
only exponential bulges embedded in otherwise normal outer disks (pentagons) are considered. 
The i? x / 4 -law bulges lie on the (7} versus absolute bulge magnitude relation traced by early-type 
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galaxies of similar luminosities, i.e. the fainter the bulge, the steeper its nuclear stellar cusp. By 
contrast, exponential bulges (extend to fainter magnitudes and) have significantly shallower stellar 
cusps for a given luminosity. It is true that these typically underly the light contributed by a 
central compact source, and are therefore subject to a larger uncertainty; however, the dichotomy 
remains even within the most conservative assumptions for the error bars. The effect is likely 
not one of "evolution" along the Hubble sequence, since we compare mostly galaxies of similar 
and intermediate, i.e. Sb to Sc, Hubble type (see Figure 5, where we plot the average nuclear 
slope inside 0.1"-0.5", (7), versus Hubble type, from RC3). In Figure 5, NGC 488 is included as 
a ii^-law bulge, while NGC 3928, ES0548G29, ESO240G12 and ES0482G17 are excluded, due 
to their uncertain morphological classification. The six squares in Figure 5 (absent in Figure 4) 
correspond to those galaxies for which no good fit to a bulge component could be derived, but for 
which the measurement of (7) could be performed (all galaxies - excluding NGC 3928 - labelled 
with "NGF" in the last column of Table 1). 

A main question is whether the difference in nuclear slopes between i? 1/,4 -law and exponential 
bulges carries important information about the structure (and formation) of these systems, or 
is rather a trivial result which simply states the difference between an inward extrapolation of 
an i? 1 / 4 -law and an exponential profile. The bulge parameters presented in papers I and II were 
obtained from fits in the range 0.5" to the last measured point: using the largest radial range 
available was found to provide an overall better fit to the data. However, performing the i? 1 / 4 -law 
and exponential fits after excluding the data inside the innermost 1" demonstrates that there is 
no ambiguity outside such radius between i? 1//4 -law and exponential bulges: the classification of 
a bulge as an or an exponential structure holds on a radial range entirely different from 

the 0.1" -0.5" used in the derivation of the average nuclear slopes (7} (see Appendix B for further 
details). This strongly supports the interpretation that esponential bulges have intrinsically 
shallower cusp slopes than i? 1 / 4 -law structures. In Figure 6 we show the L95 fits (dotted lines), 
and the single bulge or bulge plus disk fits (solid lines), for four representative objects, i.e., a single 
exponential (E482G17) or J R 1 / 4 -law (NGC488) bulge, a two-exponentials galaxy (E498G5), and 
a R l l A -\aw plus exponential object (NGC2344). When a two-components fit is required, we also 
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plot the bulge (short-dashed lines) and disk (long-dashed lines) components separately. The total 
(bulge plus disk) fits and the bulge components overlap at small radii. These examples clearly 
illustrate that while the L95 fits are aimed at representing the innermost galactic regions, the 
bulge (plus disk) fits are aimed at matching the outer profiles. 

The deprojected stellar density at 0.1" versus absolute magnitude of the bulge component 
Mfqoqw is shown in Figure 7. The symbols (as in Figure 4) represent the measurements derived 
from the deprojections of the analytical fits of Table 2 {p mo d,o.i" m Table 3). For those galaxies 
whose deprojection of the surface brightness profile gives a stellar density Pdatafi.i" which is 
significantly larger than p mo dfi.i", an arrow is plotted, which connects the value of p mo d,o.i" to the 
corresponding value of pdatafi.i"- The spherical deprojections might underestimate the density for 
large intrinsic flattening. Assuming that the exponential bulges have flattening comparable with 
that of the i? 1 / 4 -law bulges, these structures have nuclear stellar densities underlying the central 
compact sources which are at the low-end side of those of early-type spheroidals and i? x / 4 -law 
bulges. Exponential bulges might therefore break the general trend of "classical" spheroidals, i.e., 
to have a higher density the lower the luminosity. The "dichotomy" in nuclear properties between 
exponential and i? 1 / 4 -law bulges constitutes an additional piece of the puzzle concerning "how 
and when" bulges form. A discussion on the possible implications of our results for scenarios of 
bulge formation (generally supporting the idea that there is present-day "bulge-formation" in the 
centers of spirals through bar-driven infall of dissipative material; e.g., Combes et al. 1990; Hasan, 
Pfenniger & Norman 1993; Norman, Sellwood & Hasan 1996) is beyond the scope of this paper, 
and is presented in Carollo (1998). 



5. Conclusions 



In this paper we have investigated the relation between the nuclear structure of spiral galaxies 
and the physical properties of their bulges. In particular, we have (i) modelled the WFPC2 F606W 
nuclear surface brightness profiles of 41 spiral galaxies with the analytical law introduced by 
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L95 (data from papers I and II), and (ii) deconvolved the surface brightness profiles and their 
analytical fits in order to estimate the nuclear stellar densities of disk galaxies. 

Our main result is that R l l 4 -\aw bulges and exponential bulges have significantly different 
nuclear stellar cusps and densities. Specifically, i? 1 / 4 -law bulges have steep stellar cusps which 
steepen with decreasing luminosity; furthermore, their stellar cusp slopes and densities are similar 
in values to those of early-type systems of comparable luminosity. By contrast, in exponential 
bulges, the inward extrapolations underlying the light from the compact sources which sit in 
their very centers imply rather shallow stellar cusps and, very likely, relatively low nuclear stellar 
densities. 
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Appendix A. Details on the Analytical Fits 

We used the same fitting procedure described in C97a; we addess to this reference for further 
information. This was carried out in two steps: (1) The first step isolated primary from secondary 
minima of x 2 - The x 2 values were computed on a grid of points uniformly distributed on a wide 
hypercube in parameter space (with dimension equal to the number of free parameters). Once 
a minimum value was found, a new, smaller, hypercube was placed on that location, and the 



-15- 



procedure iterated. (2) The minimum value found on the hypercube was then used as starting 
point to initialize a downhill simplex minimization. We tested the procedure on simulated data, 
and verified that it recovered the initial values with high accuracy. We accepted as final the fits 
associated with the absolute minimum of x 2 . 

Deconvolutions of WF/PC data have been proven to be very reliable (e.g., L95). However, in 
our analysis, we chose not to apply any deconvolution to the post-refurbishment WFPC2 images, 
and to correct for PSF-blurring while modeling the light profiles. Therefore, the models were 
convolved with the appropriate PSFs before being compared to the data. Since pointlike sources 
with adequate S/N located near to the nuclei were not available for most of the galaxies, we 
computed the PSFs by running Tinytim (Krist 1992). Focus drifts and breathing modify the PSF 
profile and affect the flux within a 1 PC pixel radius up to 10% (and within 5 PC pixels up to 
5%; Suchkov & Casertano 1997). Therefore, the simulated PSFs obtained by construction at the 
nominal focus position are in principle of similar quality than PSFs derived from archival stars. 
Furthermore, our approach of convolving the models rather than deconvolving the data minimizes 
the effects of using a possibly non-perfect PSF. 



Appendix B. Classifying i? 1//4 -law or exponential bulges outside 1" 

In order to ensure that the classification of a bulge as an or an exponential structure 

is valid on a radial range entirely different from that used in the derivation of (7) (equal to 
0.l"-0.5"), we performed as a test the i? 1/,4 -law and exponential fits after excluding the data inside 
the innermost 1". This value is a compromise between a radius large enough to exclude the range 
where (7) is computed, and small enough to still allow the detection of the small, disk-embedded 
exponential bulges (pentagons in the figures). As an example, the results of the test are illustrated 
in Figure 8 for the same four galaxies presented in Figure 6. The solid lines represent either 
a single exponential (left panels) or a double exponential (right panels) fit; the dashed lines 
represent either a single i? x / 4 -law (left panels) or an i? x / 4 -law plus exponential (right panels) 
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fit. Two different scales are used for the abscissa for the galaxies in the left and right panels, 
consistently with the different scales of their bulge components. An offset of two magnitudes has 
been applied to ES0482G17 for plotting purposes. There are two important points to note: (i) 
the kind of profile which provides the bulge classification given in papers I and II, i.e. exponential 
or i? 1 / 4 -law, still provides a better fit to the inner galactic regions, even when the innermost 1" 
is excluded from the fits; (ii) the alternative profile with respect to the one that provides the 
classification given in papers I and II generally provides (not only a worse fit but also) physically 
meaningless best fit parameters (e.g., for ES0482G17, the i? 1//4 -law best fit of Figure 8 has an 
R e ~ 385"). We conclude that the distinction between i? 1//4 -law and exponential bulges holds in 
a radial range which excludes the one used to derive (7), and that the difference in nuclear cusp 
slopes (7} between i? 1//4 -law and exponential bulges has a physical origin. We retained the bulge 
parameters presented in papers I and II for our discussion, since those fits gave an overall better 
description of the profiles. 
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Figure 1. The surface brightness profiles, and the fits with the PSF-convolved analytical 
models (for the galaxies for which they could be derived, solid lines). 

Figure 2. The effect of the WFPC2 PSF on the deprojected density at 0.1". For ES0482G17, 
we plot the relative density error (pno-psf-Ppsf) I Ppsf versus the radius of a simulated 
compact source. Pno-psf is the density derived by deprojecting the L95 model added of the 
Plummer-modeled compact source, i.e., is the PSF-free measurement, ppsf is the density derived 
by deprojecting the profile obtained from the previous one after convolution with the WFPC2 PSF. 

Figure 3. Stellar density profiles obtained by deprojecting the observed surface brightness 
profiles (dashed lines) and their best fits (solid lines). 

Figure 4. Average nuclear slope (7) inside 0.1"-0.5" versus absolute F606W magnitude of 
the spheroidal component. Triangles are the single exponential galaxies, pentagons the double 
exponential bulges, circles the i? 1 / 4 -law bulges, and asterisks the early-type galaxies of L95 and 
F97 (small symbols) and C97a, C97b (large symbols). Conservative error bars are shown in the 
upper left corner of the diagram for the galaxies which host a central compact source (large), and 
for the remaining galaxies (small). 

Figure 5. Average nuclear slope inside 0.l"-0.5" ((7)) versus Hubble type (from RC3). 
Triangles are the single exponential galaxies, pentagons the double exponential bulges, circles the 
i? 1 / 4 -law bulges, and squares the galaxies for which no bulge-component could be fitted. 

Figure 6. The L95 fit (dotted line), and the single bulge or bulge plus disk fits (solid line), 
for four representative galaxies. E482G17 is fitted with a single exponential, NGC488 with a single 
i? 1 / 4 -law. For the two-components falaxies (E498G5, two-exponentials, and NGC2344, _R 1//4 -law 
plus exponential), we plot the bulge (short-dashed line) and disk (long-dashed line) components 
separately. The total (bulge plus disk) fits and the bulge component overlap at small radii. 
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Figure 7. Stellar density at 0.1" versus absolute F606W magnitude of the bulge component. 
Symbols are as in Figure 3. For our sample, they represent the measurements relative to the 
deprojections of the analytical fits to the surface brightness profiles (p m od,o.i") ■ For the galaxies for 
which the stellar density at 0.1" obtained from the deprojection of the observed surface brightness 
profile (pdata,o.i") ls significantly larger than p mo d,o.i"> an arrow is plotted, which connects the 
value of Pmodfi.i" to the corresponding value of Pdata,o.i" ■ 

Figure 8. Exponential and i? 1 / 4 -law fits performed as a test outside 1". The results for the 
same four galaxies of Figure 6 are presented. The solid lines are either a single exponential (left 
panels) or a double exponential (right panels) fit; the dashed lines are either a single i? 1/,4 -law (left 
panels) or an _R 1//4 -law plus exponential (right panels) fit. An offset of two magnitudes has been 
applied to ES0482G17. The kind of profile which provides the bulge classification given in papers 
I and II still provides a better fit to the inner profiles, even when the innermost 1" is excluded 
from the fits. 
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Table 1: Parameters for the 43 galaxies. Distance D and total B magnitude are from the RC3 catalogue 
(de Vaucouleurs et al. 1991). Morphological type is from the ESOLV (Lauberts & Valentijn 1989) or 
UGC (Nilson 1971) catalogues (left), or from the RC3 catalogue (right). The information provided in the 
other columns is extracted from papers I and II. £ A "Y" indicates the presence of a central compact 
source, a "N" stands for no central compact source. 1 The central source in ESO 240 G12, NGC 1483 and 
NGC 3259 is unresolved. < v ) A "Y" indicates the presence of a bar, a "N" stands for no barred component. 
A question mark indicates that a faint bar might be present, or that the entire galaxy has a flattened, 
bar-like (irregular) morphology. "ETB" indicates the presence of an early-type bar. 4 The column lists 
the total apparent magnitude in F606W and the half-light radius R e (in ") of the bulges. Close to the 
galaxy name, the superscript R 1 ^, expol or expo2 identifies galaxies whose bulge-like component has 
been fitted with a i? 1 / 4 -law or an (single or double) exponential profile. All magnitudes are corrected for 
Galactic extinction using the values published by Burstein & Heyles (1984). 
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Table 2: The best fit values for Rb, a, (3, 7, fib = —2.5 log It and (7). The column "Range" gives the 
radial interval over which the fits have been performed. When two fits are reported for the same galaxy, 
the upper values correspond to the dashed lines in Figure 1, and the lower values correspond to the solid 
lines. The average mean of the two corresponding values of (7} is used in the analysis. * indicates that 
no formal fit was performed, and the corresponding parameters optimized by visual inspection. 
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28 


N2196 


498 


900 


491 


1272 


N2344 


1156 


631 


1207 


623 


N2460 


981 


1631 


1016 


1689 








807 


1305 


N2758 


79 


193 


3 


3 


N3045 


31 


87 


6 

22 


6 

46 


IN oUu4 


318 




317 


713 


N3259 


251 


495 


107 


182 








131 


232 


N3277 


1559 


1849 


1354 


1620 


N3455 


267 


221 


20 


19 








18 


17 


N3898 


1311 


1452 


1578 


1784 


N3900 


1481 


2531 


1629 


3006 








1586 


2883 


N3928 


3900 


2199 


486 


371 


N4030 


861 


861 


151 


151 


N406 


3 


4 


2 


2 


N4260 


84 


191 


176 


272 


N4384 


6 




7 


7 


N4806 


63 


324 


1 


1 


N488 


301 


322 


328 


783 


N4980 


1 


2 


2 


2 


N5121 


1704 


2112 


1647 


2097 








2005 


2575 


N5985 


129 


485 


125 


542 


N6340 


1425 


1192 


1593 


1335 


N6384 


224 


350 


62 


71 


N7280 


1829 


2906 


1594 


3041 








1800 


3496 


N7421 


513 


859 


218 


346 








256 


412 


N7690 


1184 


1315 


1096 


1205 



Table 3: Stellar densities at 0.1" and 10 pc obtained from deprojecting the surface brightness profiles 
(subscript data) and the best fitting models of Table 2 (subscript mod). For the galaxies for which two 
best fits are available, the corresponding values for p m od,o.i" and p mo d,i0pc are reported following the same 
order as in Table 2. 
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